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Abstract 

In this talk, I discuss the general theory of neutrino oscillation experiments, putting special 
emphasis on the momentum distribution of the incoming neutrino beam. Then I discuss recent 
neutrino oscillation experiments, viz., LSND, KARMEN and CHOOZ. Experiments foreseeable 
in the near future have also been discussed at the end. 



1 Introduction 

Static properties of any particle include its mass, charge, spin, magnetic moment etc. For neutrinos, 
the charge is believed to be zero from charge conservation, and few questions are raised about this 
conclusion. There are models in which neutrinos can exhibit a small charge, but experiments have 
not paid attention to them. Similarly, the spin is known to be |, and no one thinks it would be 
interesting to challenge this value. The main unknown property of neutrinos is their masses. For 
various good reasons, the question of mass has become the obsession of neutrino physics in the last 
few decades. 

Information about neutrino mass can be obtained in a variety of laboratory experiments. We 
can divide these experiments into some distinct categories. For this, we first need to mention that in 
the standard model of electroweak interactions, neutrinos are taken to be massless particles. Thus, 
looking for neutrino mass is to go beyond the standard model. This can be done in two ways. First, 
one can study processes which are allowed by the standard model, but try to look for deviations 
in their rate from that predicted by the standard model with massless neutrinos. These are called 
kinematic tests, and these need not be accelerator or reactor experiments. For example, one can 
try to find the mass of the electron neutrino from the electron spectrum in a beta-decay process, 
which will not involve any accelerator or reactor. Thus, such tests are not the topic of this talk, 
although they will be mentioned at the end. 

A second method would be to look for processes which are forbidden in the standard model, but 
would be possible if neutrinos have masses. Here again, we want to make two subdivisions. One 
type would be processes which would require neutrino mass but not necessarily neutrino mixing. 
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Examples of such kind are neutrinoless double beta decay or neutrino magnetic moment. Again, 
these are not accelerator or reactor experiments, and so will not be discussed except at the very 
end. 

Finally, there would processes which are not allowed unless there is neutrino mixing, e.g., neu- 
trino oscillations. These use neutrino sources from reactors and accelerators, and we are going to 
describe the theory of these experiments, and some very interesting recent results. We will adopt a 
level of a person with some knowledge of weak interaction physics, but not a great deal of knowledge 
of experimental designs. This choice is not dictated by the composition of the audience. In fact, 
this reflects the level of knowledge of the speaker. 



2 Theoretical discussion of oscillation experiments 

We said earlier that neutrino oscillation is a kind of phenomena which cannot occur without neutrino 
mixing. Let us then start by explaining what neutrino mixing is. 

Consider neutrinos created in a charged current experiment, where either the initial state had a 
charge lepton i or the final state has, apart from the neutrino, a charged antilcpton i. Let us call 
this neutrino vn. For example, if one thinks of an inverse /?-dccay process where the initial state was 
e~ +]?, the final state neutrino will be called v^,. If one thinks of tt^ decay, if the charged particle 
in the final state is the neutrino would be called v^. 

This is, in fact, the sense in which so far we have used the names i/g, etc. But there is a 
catch. No one told us that the neutrino states defined this way would be the eigenstates of the 
Hamiltonian. Suppose they are not. The eigenstates will be some linear combinations of these 
states. Calling these eigenstates by i^i, V2 etc, we can write the linear superposition as 

i^e ^ / cos6l sin6' \ / z/i \ 

\v^) \ -sine COS0 A ^2 y ' ^ ■ ^ 

where ui, 1^2 etc are eigenstates. If this is the case, we would say that the neutrinos are mixed, and 
the phenomenon would be called neutrino mixing. More generally, for multiple neutrino states, we 
should write 

i^t^J^Ueai^a, (2.2) 

a 

where the subscript £ would run over the states Ue, t'^, i^r etc, and z/q's would stand for the eigen- 
states. However, for the sake of simplicity, we will work in a two-level system. 
Suppose now we have created an initial beam 

|i/(0)) = li/a) — cosaji^e) -|- sin a 

= cos(6l + a) + sm(e + a) \u2) . (2.3) 

After this beam travels for a time t, the state would evolve to 



\u{t)) = e-'"' = e-'^'' cos{9 + a) + e''^^' sm{9 + a) \u2) 

'cos(e + a) + e-'^^^ sm(e + a) \u2)] , (2.4) 
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where 

AE^E^-E,^^, (2.5) 

where K is the magnitude of the 3-momentuin. In writing the last step, we assumed that the 
neutrinos are ultra-relativistic, so that the energy-momentum relation can be approximated by 

ui^ 

E^^flO^^^K^ — , (2.6) 

neglecting higher powers of mass. 

Now suppose we are trying to find the state vp, which is defined similar to except the angle 
a replaced by the angle (3. The probability for this would be 

\{vii\v{t))\^ = cos(^ + a)cos(^ + /?) + e-^^'^*sin(^ + a)sin(^ + /?) ^ 



1 + cos(2^ + 2a) cos{2e + 2(5) + sin(2^ + 2a) sm{2e + 2(5) cos ( -i^x 



1 /Am^ \ 

= - [1 + cos(2a - 2(5)] - sin(2^ + 2a) sin(2^ + 2(5) sin^ ( -^x j . (2.7) 

replacing t by a; since the neutrinos are ultra-relativistic anyway. 

In particular, if we consider a — (5, i.e., we try to find the same state that was created ect t — 0, 
we obtain the "survival probability" of that state to be 

Psurv = 1 - sin2(2^ + 2a) sin^ (^^^x^ . (2.8) 

On the other hand, if is orthogonal to the original state |i/q,), we obtain is the "conversion 
probability", which is 

(Am^ \ 

If we consider only the pure flavor states, i.e., or i/^ for the initial and final states, a is either 
or 7r/2. In both cases, we obtain 



/ Am 
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-Psurv = 1-sin 26' sin y-j^^ j > 

(Am^ \ 
l^^j ■ ^^'^^^ 

To connect this formula to the results of real experiments, we first need to realize that in a given 
experiment, the neutrinos can never come with a well-defined value of K due to the uncertainty 
relation. Thus, the conversion probability obtained in a real experiment should be given by 

dK^K) sin^ [^x \ , (2.11) 

where ^{K) is the spectrum of the incoming beam, normalized by 

J dK ^{K) = 1 . (2.12) 



— 1 1 1 1 — I — I — r- 



— 1 1 1 1 1 — I — I — r 



P •> P 




Figure 1: Implication of a hypothetical neutrino oscillation experiment with a gaussian momentum 
distribution, with ax/iK) = 0.2, which obtains Pconv = 0.01, or, equivalently, Psurv = 0.99. The 
dashed line represents the approximate form of the equi-probability contour given in Eq. (|2.23| ). 



Suppose now in an experiment one obtains Pconv = Po- What will it mean? ^{K) should be 
known from the design of the experiment and the nature of the source. One is then left with two 
variables, sin^ 2^ and Am^. For reasons that would be clear as we proceed, it is useful to use the 
dimensionless quantity 



Am? 
Ta = —r^TT X 



4{K) 

instead of Am^ directly. The conversion probability can then be written as 



sin^ 29 J dK $(7^) sin^ Ir^ 



K 



(2.13) 



(2.14) 



Thus, the result of this experiment would pick up an allowed contour in a plot of ta vs sin^ 29. The 
results for Pq = 0.01 are shown in Fig. using a gaussian distribution of momenta: 



(TkV2tt 



exp — 



2al 



(2.15) 



where we have taken ax = 0.2{K) for the sake of illustration. 

Let us discuss the asymptotic features of this figure. If ta is large, it is useful to rewrite Eq. ( p.l4|) 
in the form 



Pconv (x) = ^ sin'' 29 J dK^K) 



1 — cos 2rA 



K 



(2.16) 



The argument of the cosine function oscillates violently as a function of momentum in the region 
where $(i^) is appreciably different from zero. Thus, the integral of this part vanishes. The other 
part gives 



sin^ 29 = 2Po for ta > 1 , (2.17) 
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using Eq. ( p. 12 ). This is the vertical part of the hne at the top of Fig. |l|. 



On the other hand, if is small, we can replace the sine function by its argument. Eq. ( p.l4| ) 
now takes the form 



Po = rlsin^ 29 {K)^ [ dK <I>{K) — 

J K 



This can be rewritten as 



where 



rlsin' 29 {Kf l^-^) . (2.18) 



ri sin^ 29 = for ta < 1 , (2.19) 



p 

The denominator of the right side of this equation is determined by the momentum distribution. 
Therefore, Pq is known from the results of the experiment. Eq. ( |2.19| ) now says two things. First, 
the slope of the equi-probability contour in the log-log plot of ta vs sin^ 29 should be — |, as seen in 
the lower part of Fig. |l|. Second, when sin^ 29 = 1, the value of ta for the equi-probability contour 
will be given by 

rA = for sin^ 29 = 1. (2.21) 

This, in fact, is the lowest value of ta probed by the experiment. 

Since small ta means more sensitive information about masses, let us spend a little bit of time 
on this formula. Consider what happens if we have a momentum distribution with very little spread. 
In this case, we can write {l/K"^) ^ 1/ (K)"^, so that the denominator of Eq. (|2.2CI|) would be unity. 



In this case, Pq ~ Pq, so that the lowest value of ta probed by the experiment is really This 

is the case with the curve of Fig. [l|, for which a small momentum spread was assumed. However, 
in general, Pq 7^ Pq. Rather, 

Po < Po (2.22) 

since in any distribution of a non-negative variable, (K)"^ {l/K"^) > 1. This is because larger values 
of the variable will be preferentially picked up by the average, whereas smaller values will be favored 
in the evaluation of {l/K"^). Thus, the spread of momentum of the initial beam is not a hassle that 
one has to cope with in these experiments. On the contrary, it is helpful to have a good spread, 
which allows one to probe lower and lower values of ta- 

Notice that although we used a specific form of the momentum distribution for plotting Fig. |T|, 
we have never used it in discussing the asymptotic behavior of the plot. Thus, we can say that 
no matter what the momentum distribution is, we would obtain the asymptotic forms given in 
Eq. (|2.17|) and Eq. ( p.l9| ). Of course, the value of Pq will depend on the distribution, and so will 
the wiggles in the plot which appear in the region around ta ~ 1, between the two asymptotes. If 
we do not care about these wiggles, the rest of the equi-probability line can be represented by the 
analytical formula 

sin^ 29 = 2Pq + ( ^ - 2Po ) exp(-ri) . (2.23) 




This has also been shown in Fig. |T] with a dashed hne, taking Pq = Pq in view of the low standard 
deviation used for the plot. 

Of course in real experiments one does not obtain the conversion or the survival probability 
without any error bar. If the experiment is a null experiment, i.e., if they find that the conversion 
probability, e.g., is less than some value, this excludes the upper right portion of the parameter 
space. On the other hand, if the experiment obtains a positive signal, we need to draw two equi- 
probability lines corresponding to the smallest and the largest probabilities allowed, and the region 
between two lines will be allowed. 



3 Different classes of oscillation experiments 

We now classify the existing experiments into various categories for the clarity of presentation. The 
classifications will be done on various grounds. 

Values of (K) : The value of the energy of the neutrinos is determined by the type of sources of 
the neutrinos. For reactor neutrinos, the neutrinos are produced as by-products of nuclear reactions, 
typically fission reactions, and so the energies are of the order of a few MeV. The energies are not 
high enough to produce the charged muon or tau, so no i/^ or Uj. is available in this source. Moreover, 
since fission reactions produce nuclei of lower atomic number in which the neutron to proton ratio 
is usually lower, at the nucleon level these reactions convert neutrons to protons via the reaction 
n ^ p + e~ + Ue- Hence the neutrino source is z/g for this kind of experiments. Most of the early 
experiments are of this type. Recently, the CHOOZ experiment have given their first results. We 
will discuss the results later. 

Another kind of experiments use medium energy neutrinos, which come in the range of a few 
tens of Me Vs. The sources here are typically pions created by beam dump of a proton beam. At 
present, there are two oscillation experiments which use such neutrinos. They are KARMEN and 
LSND. 

High energy accelerator neutrinos, with energies in the range of GeVs, are also beginning to be 
used. Many such experiments are also in the planning or constructing stages, and we will mention 
some. 



Values of x : Early experiments on neutrino oscillation used path lengths of the order of at 
most a few hundred meters. As indicated above, any given experiment, with a given accuracy of 
detection, will probe the parameter space only upto a certain value of ta which is roughly equal to 
\fP^- So, to probe lower values of Am^, one needs large x. Thus, some long baseline experiments 
are also being planned. 

To obtain a feeling for the magnitudes, it is useful to use the form 

— (^)(t^)""(i^)- 

Thus, for example, if one wants to probe Am^ values at the level of 10~^ eV^ with an experiment 
where Pq ~ 10~^, one needs x/ (K) ~ lOm/MeV. This must be taken only as a benchmark value, 
rather than a firm estimates, because we used \/Po as the minimum value of ta that can be probed. 
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As we remarked earlier, this is true only if the standard deviation of the beam is very small compared 
to the mean value. In real experiments, this is not often the case. 

Appearance and disappearance experiments : The experiments either look for the survival 
probability of a certain kind of neutrinos, or the conversion probability of a particular flavor to 
another particular flavor. The first kind of experiments are called disappearance type, and the 
second kind appearance type of experiments. 

If there were only two kinds of neutrinos, the two approaches would provide the same informa- 
tion. But there are at least three kinds of neutrinos, so this is not true. There are advantages and 
disadvantages of both types. The appearance experiments are easier in the sense that it is easier to 
detect a flavor which was not present in the original beam — even one event would signal a non- 
trivial signal. In the disappearance experiments, on the other hand, one has to check meticulously 
whether the flux of the original beam has decreased over the distance. 

However, appearance experiments can check only a specific channel. For example, suppose we 
try to do an experiment to see whether z/g oscillates to z/^, and we find a null signal to the accuracy of 
the experiment. So z/g does not oscillate to u^, but who knows, maybe it still oscillates substantially 
to Ur, or some other unknown variety of neutrinos. If on the other hand we try a disappearance 
experiment of an original z/g beam, we would get a non-trivial result no matter what z/g oscillates 
to. But it is also true that in this experiment, we will not know what z/g really oscillates to. 

4 Recent results of oscillation experiments 

We now present the results of some recent experiments, and prospects for some upcoming ones. 
LSND 

The Liquid Scintillator Neutrino Detector group, working at the Los Alamos Meson Physics Facility 
(LAMPF), have consistently reported some non-null results over the last few years. If confirmed, 
these will constitute the first positive evidences for neutrino oscillations. 

For the neutrino source, they produce pions by dumping 800 MeV protons on a target. Due to 
charge conservation, vr"*" is produced predominantly. Then vr"*" decays to /i"*" + z/^, and subsequently 
/i"*" decays to e+ + z/g -|- z/^. Thus, they have z/g, z/^, and z/^ in equal amounts. 

The first experiment of this group sought for z/^ Vf. oscillations, with the coming from 
7r+ decays at rest. Recently, they have published [Q] their results for z/g oscillations, using 

decays of pions in flight in order to obtain a broad momentum distribution. We show both results 
in Fig. |. 

A few things can be noticed from their plots. First, since they used in-flight pions for the 
—>■ oscillations, the spectrum was very broad, and consequently the wiggles have been almost 
smoothed out. Second, note that CPT invariance tells us that the mixing angles and Am^ for both 
these cases should be the same. Indeed, the regions allowed by the two have a lot of intersection, 
which serves as a good cross check for the results. 
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Figure 2: Results of the two LSND experiments. The dashed lines indicate the allowed region from the 
t'/i Ve experiment, and the solid lines the same for the — > v^. experiment. From Ref. 0. 



Table 1: Details of the LSND and KARMEN experiments. 





LSND 


KARMEN 


Source 


Beam- dump of 800 MeV pro- 
tons from linear accelerator 


Beam-dump of 800 MeV pro- 
tons on a Ta-D20 target 


Production 
Mechanism 


7r+ ^ + z/^, and subsequently fi'^ + Ue + i'^ 


Produced flavors 


z/e, z/^, z/^ (in equal amounts) 


Energies 


z/^ : 30 MeV (mono-energetic) 

z/g and z/^j : to 53 MeV (different dispersions) 


Po for v^, Ue 
Po for Uf, Ue 


(2.6 ±1.0 ±0.5) X 10^3 
(3.ltl:J±0.5) X 10"3 


< 2.0 X 10-2 

< 4.3 X 10-3 
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Figure 3: Results of the KARMEN experiment. Also shown are the expectations for an upgraded 
experiment, and the allowed regions from some earlier experiments, including the first LSND experiment. 
From Ref. |[]. 



KARMEN 

This is the KArlsruhe Rutherford Medium Energy Neutrino experiment 0. The experiment is very 
similar to the LSND one, in the sense that they produce neutrinos by beam dump as well, although 
the dump material is not the same. The two experiments, LSND and KARMEN, are summarized 
in Table |l[ 

The KARMEN group also, like the LSND group, have performed two kinds of appearance 
experiments, z/^ —>■ v^. and z/g. The second one is more sensitive, since the source has no 

Vf.. They do not find any indication of oscillation in either experiment. They find only upper 
limits of conversion probabilities. The results are shown in Fig. |^, where the results of some earlier 
experiments are also shown, including those from a reactor experiment from the Bugey reactor in 
France, and a Brookhaven experiment BNL E776. 

In the figures, one can see the region allowed by the LSND experiment of z/^ — > z/g. As is seen, 
there is no disagreement between the two results at the moment. However, the KARMEN group 
is working on an upgrade. With this, they will probe much lower values, which are shown in the 
figure as well. With this upgrade, they should see positive results if the LSND results are correct. 
And if they don't see anything, that will be a contradiction with the LSND results. 



CHOOZ 

Very recently, results have come from a reactor experiment in France |Q. The experiment is done 
with (i^) ^ 3 MeV and a; ~ 1 km. Since the source of neutrinos is a reactor, the primary beam 
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Figure 4: Results of the CHOOZ experiment. Also shown are the results of some earlier reactor experi- 
ments, and the region favored by atmospheric neutrino anomaly data from Kamiokande. From Ref. [ffl]. 



consists of Ve- And the experiment is of disappearance type, looking for v^. ~^ i-'X) where X can be 
any flavor, including unknown ones. They obtain 

l-Po = 0.98 ±0.04 ±0.04. (4.1) 

This is consistent with no oscillations at all. At 90% CL, they use therefore 

Po<0.09. (4.2) 

The important thing is that, even with a rather high value for Pq, they can probe quite small 
values of Aim? since {K) is small and the path length x is large. Their results are shown in Fig. ^ 
where results of some earlier reactor experiments have also been shown. 

In addition, we have also shown the region of parameter space that can explain the atmospheric 
neutrino anomaly measured by the Kamiokande collaboration, provided the anomaly is caused by 
oscillations. The CHOOZ result has ruled out the region. This of course does not mean 
that either CHOOZ or Kamiokande is wrong. It means that, if the atmospheric neutrino anomaly 
has to be explained by neutrino oscillations, it must be because the z/^'s are being lost by oscillating 
to or to some unknown species. 



5 Outlook for future experiments 

Let us reiterate the statement made earlier, that in fact an oscillation experiment measures directly 
the parameter ta. So, with a given sensitivity of an experiment (i.e., given Pq)? one can probe 
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smaller and smaller values of Am^ if one can construct experiments with larger values of the path 
length X and smaller values of the average energy of the neutrino beam. It would thus seem that 
if we can perform reactor experiments with path lengths as large as possible, we would obtain the 
best results. 

There is however a limitation to this goal. With reactor neutrinos, first of all, one can do dis- 
appearance experiments only, and measuring the survival probability is difficult. Secondly, reactor 
beams are not very well directed. So, if we try to go to large path lengths, we also lose considerably 
in the sensitivity. One can already see it in the results of the CHOOZ experiment, which uses a 
path length of 1 km which is considerably higher than those of earlier reactor experiments. But 
they paid the price by ending up with a Pq which is considerably higher compared to these earlier 
experiments. One can see this from the fact that in Fig. ^, the vertical part of the line for large Am^ 
occurs at much higher value of sin^ 29 for the CHOOZ experiment compared to the earlier ones. 
With reactor experiments, it is hard to go much beyond in terms of path length. Nevertheless, the 
Kamiokande group is contemplating an experiment with neutrinos from nearby reactors, where the 
path lengths could well be larger than a kilometer. 

Accelerator beams are however well focused and one can hope to go to hundreds of kilometers 
of path length with them. Various long baseline experiments are in the planning or construction 
stage at the moment with accelerator neutrinos. We give a brief summary of them 

K2K : The distance between KEK laboratory and the Kamiokande experiment in Japan is 
about 235 km. There is a plan of sending a neutrino beam from KEK to Kamiokande, with an 
average energy of 1 GeV. At Kamiokande, they will be detected. The beam-line should be finished 
by 1998, and data should be expected by 1999. 

Fermilab— Soudan : The distance here is 735 km. The MINOS experiment will be installed 
in the Soudan mines. The project may start at the beginning of the next century. 

CERN— Gran-Sasso : The distance is 732 km. Neutrino beam from CERN can be detected 
at Gran Sasso. Various proposals exists for detection of the neutrinos at Gran Sasso. 

6 Conclusion 

Since I have started with an "Introduction", I have to end with a "Conclusion". There are of 
course some obvious concluding remarks. With more experiments, we will know things better. The 
important point is that at the moment, there seems to be a lot of activity in the field. With the 
KARMEN upgrade, the LSND results will be checked, and we are looking forward to it. The long 
baseline experiments will usher a new era in this field, although the results of these experiments are 
not expected until a few years from now. 

One of the things that I want to emphasize here is that the neutrino oscillation experiments at 
present provide the most sensitive information about neutrino mass from all terrestrial experiments. 
Let us compare them briefiy with other kinds of experiments described in the Introduction. The 
kinematic tests provide only upper bounds of neutrino masses, which is of the order of a few eV for 
Ue, 170 keV for i/^, and 24MeV for u^. 
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Neutrinoless double beta decay process is not allowed unless lepton number symmetry is violated. 
The amplitude of this process depends on the factor 

meS = ^U^a^a, (6.1) 

a 

where is the mass of the a-th neutrino eigenstate. In absence of neutrino mixing, the quantity 
rries becomes equal to the mass of the z/g. The most stringent bounds come from the Heidelberg- 
Moscow experiment 0. They use 19.2 kg of ''^Ge and search for the process 

''Ge^ ''Se + 2e-, (6.2) 

and set the following lower limit on the lifetime: 

Ti/2 > 7.4 X 10^^ yr. (6.3) 

It is not straight forward to deduce the upper bound on nies from this raw experimental number. 
The reason is that the calculation of this rate involves nuclear matrix elements, the evaluation of 
which cannot be done in a manner that pleases everyone. Using various evaluation of the matrix 
elements, one obtains the upper bound on ranging between 0.56 eV to 1.76 eV. One can say 
that the upper limit is of the order of an eV. 

Neutrino oscillation experiments, on the other hand, cannot probe the mass values directly. 
They can measure the mass squared differences. If, inspired by the mass patterns of the charged 
fermions, one assumes a hierarchy in the neutrino masses, the lowest mass value probed would be 
the square root of the lowest Am^ probed. With this assumption, we see the neutrino oscillation 
experiments are already probing mass values which are lower than an eV. With upgrades and 
improvements, they can go to even lower values. Fortunately, these upgrades and improvements are 
not day dreaming, they are things to expect in the near future. 
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